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ABSTRACT 

We discuss a possible mechanism for heating the solar atmosphere by the ensemble of ther- 
mal waves, generated by the photospheric dynamo and propagating upwards with increasing 
magnitudes. These waves are self-sustained and amplified due to the specific dependence of the 
efficiency of heat release by Ohmic dissipation on the ratio of the collisional to gyro- frequencies, 
which in its turn is determined by the temperature profile formed in the wave. In the case of 
sufficiently strong driving, such a mechanism can increase the plasma temperature by a few times, 
i.e. it may be responsible for heating the chromosphere and the base of the transition region. 

Subject headings: Magnetohydrodynamics — Sun: chromosphere — Sun: transition region 



1. INTRODUCTION 

Seeking for the mechanism of heating the so- 
lar atmosphere with height from a few thousand 
to a million Kelvin is one of long-standing prob- 
lems in astrophysics. The approaches proposed 
by now, roughly speaking, can be separated into 
two groups (e.g. review s bv I Walsh fc Irelandll2003t 
lErdelvi fc Ballail 120071 and references therein). 
The first group deals with a generation of some 
kinds of (magneto-) hydrodynamic waves or pulses 
in the base of the solar atmosphere and their sub- 
sequent propagation and dissipation in the upper 
layers. The mechanisms of the second group as- 
sume that the heating is due to the ensemble of 
small-scale flare-like events, caused by the recon- 
nection processes in the specific magnetic field 
configurations. It is commonly believed now that 
no unique mechanism can provide the entire heat- 
ing of the solar atmosphere, and a few of them are 
acting simultaneously in the Sun. 



Besides, as was proposed bv I Aschwanden et al 



(2007), the inverted (increasing with height) tem- 
perature profile might be formed "dynamically" — 



due to the fluxes of plasma heated in the lower lay- 
ers and propagating upwards (which is often called 
evaporation). This is, in fact, the third kind of the 
heating mechanisms. 

The aim of the present work is to consider one 
more mechanism, which, strictly speaking, belongs 
to none of the three above-mentioned groups but 
inherits some features from all of them. Like mech- 
anisms from group II, it assumes an ensemble of 
the small-scale flare-like events; but such events 
require no special field topology for their occur- 
rence (as in the case of reconnection) and are more 
similar to the specific type of wave pulses, as in 
the mechanisms of group I. At last, our scenario 
reminds the Aschwanden's approach (group III), 
because the heating process propagates dynami- 
cally from the lower to upper layers. 

2. FORMULATION OF THE MODEL 

The dynamo processes in weakly-ionized plas- 
mas, considered in this section, closely remind 
the so-called S^-dynamo, well k nown in the iono- 



spheric physics (e.g. review by iRichmo nd 1989, 
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Fig. 1. — Sketch of the dynamo processes in the partially-ionized plasma (motion of the neutrals, electrons, 
and ions are colored green, blue, and red, respectively, in the online version of this figure) and a schematic 
structure of the associated thermal wave, propagating upwards (colored magenta in the online version). 



and references therein). Similar processes in the 
solar physics are called photospheric dynamo, 
and they were employed in a number of pa- 
pers to describe generation of field-aligned cur- 
rents and perturbations of the magnetic fields 
(jHenoux fe Somovl 119871 ). The associated ther- 
mal effects were usually neglected; a n exception is 
the old work bv lSen fc Whitel (|l972h . 



2.1. Qualitative Description 

Let us consider a column of partially-ionized 
plasma stratified by the gravitational field g and 
permeated by the magnetic field B, which for sim- 
plicity is taken to be vertical and constant. Be- 
sides, we assume that this plasma is driven radially 
in the horizontal plane by the motion of a neutral 
component v„ (which can be associated, for ex- 
ample, with convective flows at the photospheric 
level); see Figure [TJ 

Due to the gravitational stratification, a neutral 
gas density sharply decreases with height. Respec- 
tively, the free-path length and the time of inter- 
particle collisions sharply increase, while the col- 
lisional frequency decreases. Then, at the bottom 
of the column, where the frequencies of collisions 
of both electrons and ions with neutrals, v en and 
fin respectively, are much greater than their gy- 
rofrcquencies uj e and u>i, the charged particles are 
completely dragged by the neutrals and move with 
them; so that no local electric fields and currents 
are generated. 

In the opposite case, at the top of the col- 
umn, where the collisional frequencies are much 
less than the gyrofrequencies, the electrons and 



ions are unaffected by the radial motion of the 
neutral component and, therefore, will remain at 
rest. No electric fields and currents are produced 
again. 

The most nontrivial situation arises at the in- 



termediate heights, where Vy, 



• uji and v en < uj„ 



Under these circumstances, the ions are dragged 
by neutrals, while the electrons remain approxi- 
mately at rest. Due to the separation of the elec- 
trons and ions, a radial electric field E should de- 
velop, resulting in its turn in the azimuthal drift 
of the charged particles. However, while the elec- 
trons will experience such a drift almost without 
resistance, the ions will be decelerated by the col- 
lisions with neutrals. Therefore, the azimuthal 
(Hall) electric currents j are generated, as shown 
in Figure [TJ 

The above-mentioned electric fields and cur- 
rents are concentrated in a quite narrow height 
range, where z^„ ^w^; and they should lead, firstly, 
to the magnetic-field pe rturbations (considered by 



Henoux fc SomovMl987l ) and, secondly, to the dis- 



sipative effects and heat release, which are just the 
subject of our study. (The dissipation exists due 
to the azimuthal electric field appearing in the co- 
ordinate frame co- moving with plasma.) 

Once a heat release has occurred at the height 
where i>i n » (we shall call it z e s), this spot is 
no longer efficient for the subsequent heat release, 
because the plasma is heated and the thermal ve- 
locity of its particles Uj as well as the collisional fre- 
quency Vi n = Vi j Xi n increased (under assumption 
of time-independent neutral density n n , constant 
ion-neutral cross-section oi n and, consequently, 
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the unchanged mean free path Aj„ = l/(crj n n n )). 
Therefore, the condition Vi n K,uji will no longer be 
satisfied. 

On the other hand, due to the heat conduc- 
tion in space, the ion thermal velocity u, will in- 
crease also at a greater height (where the mean 
free path X in is larger due to the initial gas strat- 
ification) and, therefore, the spot of effective heat 
release, V{ n = Uj/Aj„ rs u>i, will move there. As a 
result, we should expect formation of the specific 
self-sustained thermal wave propagating upwards, 
as depicted in the right-hand panel of Figure Q] 
Moreover, the wave amplitude will increase with 
height, just because of decreasing a heat capacity 
per unit volume in the stratified gas. 

It should be also mentioned that heat release in 
the dynamo-region may be further facilitated by 
the o nset of two-stream instability (pen fc White! 
1972L and references therein) ; but we shall not con- 
sider here in more detail the corresponding micro- 
physical mechanisms. 

Finally, let us discuss more carefully the scope 
of applicability of our scenario. In fact, the ion's 
collisional frequency is composed of two parts: 
Vi = fin + Va, where z^ n and i>n are the fre- 
quencies of ion-neutral an ion-ion collisions. It 
was assumed everywhere above that the second 
term can be ignored. On the other hand, the de- 
gree of ionization increases with height; so that 
the ion-ion collisions will inevitably become dom- 
inant starting from some altitude. In this case, 



\ H = l/(<Tum) ~ vf/rii 



since a;. 



So, 



the ion's collisional frequency vu— Vi/Xu ~ rii/vf 
decreases with temperature, as distinct from the 
case of ion- neutral collisions. However, this fact 
is unrelated to the criterion of the dynamo-layer 
development, Vin^^h because collisions between 
the charged particles of the same kind cannot af- 
fect the generation of electric currents. The dom- 
inant role of ion-ion collisions can change only the 
coefficient of heat conductivity. 

Therefore, the scope of applicability of the pre- 
sented model is somewhat wider than seems at the 
first sight. At the same time, the plasma must be 
weakly ionized in the entire height range under 
consideration, since otherwise it would be mean- 
ingless to speak about dragging the charged par- 
ticles by the neutral gas. 



2.2. Basic Equations 

In the simplest one-dimensional approximation, 
the process of heat transfer along the magnetic 
flux tube can be described by the equation: 
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where T is the temperature of the heavy particles, 
i.e. neutrals and ions (the electron temperature is 
not of interest here), t is the time, z is the ver- 
tical coordinate, n is the number density of the 
heavy particles, x ~ (1/3) Aw is the temperature 
conductivity, c = const is the heat capacity per a 
heavy particle, and q is the volume density of the 
heat release due to Ohmic dissipation. (The ba- 
sic parameters of neutrals and ions are assumed to 
be approximately the same and, therefore, written 
without additional subscripts.) 

Leaving aside a microscopic theory of the 
Ohmic heating, let us use a phenomenological ap- 
proximation of the heat-release profile, based on 
its qualitative behavior discussed above. Namely, 
we take Taylor expansion in the vicinity of the 
heat-release maximum, achieved at Vi n ~ uji. 



q{z,T) = q < 1- a 1 - 



v in {z,T) 



, (2) 



where qo is the amplitude of the heat release, and 
a is the characteristic width of the profile. 

Next, we shall use for simplicity the time- 
independent exponential (i.e., actually, isother- 
mal) height profile of the gas density: 



n(z) — hq exp[— z/H] 



(3) 



where H = T /mg = const is the height scale. (In 
other words, it is assumed that the initial gas pro- 
file is not changed during the heat propagation; of 
course, redistribution of the gas density should be 
taken into account in a more accurate treatment.) 

The thermal velocity is evidently related to the 
temperature as 



''0 



(4) 



So, assuming that the cross section a is constant 
(which is a reasonable approximation for the colli- 
sions of ions with neutrals and neutrals with each 
other), we get: 



v in = anv y/T/T = v in0 (n/n )^T/T . (5) 
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Besides, it is convenient to choose the origin of 
z-axis (z = 0) in the point where v% n Q=u)i. 

At last, substituting expressions ([2j)-([5]) into 
equation ((T|), we obtain the equation of heat trans- 
fer in the following form: 
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(6) 



where the dimensionless quantities (marked by as- 
terisks) were introduced as 



T*=T/T , t* = t/r . 
H* = H/\ , q% = 



z*=z/\ , (7) 
qoTo 



(8) 



cn T ' 

To is the initial temperature of the gas, tq= v~^ , 
and Ao = Aj„o ■ 

3. RESULTS OF THE COMPUTATION 

Expression (|6]) represents a quite specific 
type of the nonlinear parabolic partial differ- 
ential equation. A convenient tool for solving 
just this type of problems was implemented in 
the NAG nu merical library — this is the subrou- 



tine D03PCF (jNumerical Algorithms Group||2001 
Berzind[l990h 



We restrict our consideration here by the sim- 
plest numerical solution, just to illustrate how the 
proposed mechanism works. (A detailed treat- 
ment of more realistic situations will be published 
elsewhere.) So, let us take the initial condition in 
the form: 

T* = 1 at t* = (9) 

(i.e. the entire gas had initially the same temper- 
ature To). 

The boundary conditions will be specified, 
somewhat arbitrarily, as vanishing heat fluxes at 
the top and bottom of the computational region: 



dT* 

dz* 



at 



and 



(10) 



(However, we shall discuss below only behaviour 
of the solution far from the boundaries, which is 
insensitive to the boundary conditions.) 

At last, let us take, for example, the following 
parameters of the problem: H* = 10, q$ = 0.1, 




Fig. 2. — Computed height profiles of the temper- 
ature at the successive instants of time (colored 
red in the online version of the figure). The star 
(colored blue) shows the spot of the initial heat re- 
lease, subsequently moving to greater heights; and 
the dashed (grey) curve represents a hypothetical 
temperature profile that would be formed in the 
case of a static heat source. 

a = 1, z* d = -100, and z* = 100. The obtained 
numerical solution is presented in Figure [21 where 
height profiles of the temperature are drawn for 
the successive instants of time with increments 
At* = 10. 

As is seen, a propagation of the self-sustained 
thermal wave really results in the formation 
of temperature profiles monotonically increasing 
with height above the spot of the initial heat re- 
lease at z = 0. This is essentially related to the 
dynamical nature of the process; otherwise, at 
the fixed heating source, the temperature pro- 
file would be strongly asymmetric (because of the 
substantial height dependence of the heat conduc- 
tivity), but its maximum would always take place 
at the fixed site of the heat release. 

Next, the series of profiles in Figure [2] monoton- 
ically increases with time, without any clear evi- 
dence for saturation. This is just because our sim- 
plified model does not take into account the radia- 
tive loss as well as exhausting efficiency of the dy- 
namo processes in the course of time. It would be 
desirable, of course, to include into the future nu- 
merical simulations a realistic function of radiative 
loss (which should be especially important at the 
enhanced degree of ionization, associated with the 
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increased plasma temperature). However, the role 
of radiative cooling should not be overestimated 
in advance: when the plasma temperature is well 
below the ionization thresholds of the plasma ele- 
ments, then the establishment of higher ionization 
density requires a very large number of interparti- 
cle collisions and, therefore, may take a consider- 
able time. So, if the discussed thermal waves are 
impulsive phenomena, then there may be just in- 
sufficient time for the increase in ionization, and 
the radiative losses will be not so large. 

4. DISCUSSION AND CONCLUSIONS 

As follows from the above consideration, the 
presented model may be a mechanism for heating 
the lower part of the solar atmosphere, i.e. for- 
mation of the temperature profile increasing with 
height. At the same time, it should be borne in 
mind that it works only in weakly-ionized plasma. 
To provide plasma heating up to the coronal tem- 
peratures, it is necessary to invoke other heat- 
release mechanisms, commonly based on the mag- 
net ic reconnection (for t he review of recent re sults, 
see 



Shibata et al~ll2007t lAulanier et al.ll2007l ) 



Let us mention also that some qualitative ideas 
why heating of the solar plasma should occur just 
at the spot of approximate equality between the 
collisional and gyro- frequencies we re put forwar d 
in our paper about a decade ago (jDumin 
Unfortunately, it was not clear from that work 
how is it possible to get the temperature profile 
monotonically increasing with height rather than 
peaking at a fixed altitude of the maximum heat 
release. The mechanism of the propagating ther- 
mal wave, described in the present paper, gives a 
possible answer to that question. 
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